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ABSTRACT 

We review the enveloping algebra of the 10 dimensional chiral sigma matrices. To 
facilitate the computation of the product of several chiral sigma matrices we have de- 
veloped a symbolic program. Using this program one can reduce the multiplication of 
the sigma matrices down to linear combinations of irreducilbe elements. We are able to 
quickly derive several identities that are not restricted to traces. A copy of the program 
written in the Mathematica language is provided for the community. 
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I. SPINORS AND SIGMA MATRICES IN 10 D 



Several problems in 10 D Physics require the ability to reduce products of the 10 D chiral sigma 
matrices into a linear combination of the irreducible matrices [1-4]. In general such calculations 
can quickly become intractable. For some problems, many approximation schemes are applied such 
as linearization and perturbation theory. However when one is discussing properties of algebras, 
such as the closing of the Jacobi identity or commutation relations, exact results are mandatory 
[5]. Computations such as computing the equations of motion in some supergravity theories, could 
literally take several weeks. These types of calculations could be reduced to a few hours work 
using today's powerful symbolic manipulators, such as Mathematica and Maple, in conjuction with 
high speed computers. In this note wo provide the community with an example of such a program 
along with several identities that we have derived using this tool. We begin by outlining the formal 
construction of the enveloping algebra of the 10 dimensional chiral sigma matrices. 

It is well known that in D-dimensional space-time Dirac spinor has 2"? complex components if D 
is even and in that case all the representations of the Dirac algebra are equivalent, while in odd 
dimensional space-time the number of components is with two inequivalent representations. 

This result is valid for any signature of the flat space-time metric. However, this kind of spinor is 
not, for any D, the smallest irreducible representation. When D is even we can always impose the 
Weyl condition which separates the spinor into its left and right handed components which transform 
independently under the Lorentz group. Each of these Weyl spinors thus has only half degrees of 
freedom of the original Dirac spinor. One more condition can be imposed to reduce further the 
number of independent components. This is so called Majorana condition which relates the charge- 
conjugate Dirac spinor to the original one. If after the charge conjugation the spinor is not changed, 
the number of independent components is reduced again by a factor of two. In that case the left 
and right Weyl spinors are related. The resulting spinor is necessarily real. This type of reduction 
can be done in 2 (mod 8) dimensional space-time. It then follows that for D = 10 we can impose 
both conditions and have for the basic spinor a 16-dimensional real object. Consequently, instead of 
using the Dirac gamma matrices we use the Pauli matrices, of dimension 16 x 16 with real elements. 
We can then use small Greek index to denote 16 component left-handed Majorana- Weyl spinor and 
dotted index for a right-handed one: 

{rr = ir), = (x"). (i) 

The indices can be raised and lowered by the use of the charge conjugation matrix which has mixed 
indices: C^^ and C"^, 

% = rC^p^ Xa = X^C^p, (2) 
C^^C''^ = 4, C'^C^^ = 5^. (3) 

As we can see, raising and lowering the indices also changes the handedness of a spinor. The Dirac 
algebra (of the Pauli matrices) is defined through the usual anticommutation relation: 

The sigma matrices can be regarded as bispinors. There are three types of them: purely left-handed, 
purely right-handed and mixed ones. The purely left-handed ones are 
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(5) 



where lower case Latin indices are vector indices denoting space-time directions. Purely right-handed 
bispinors have only dotted indices, but due to the existence of the charge conjugation matrix we 
have 



The mixed bispinors are 



They are, of course, related to 



{a2-)°'^ = C""C'^'^{a2-)^g. (6) 



Sf, (^^)/, (8) 



The sigma matrices with more vector indices are defined through the following multiplication table: 

i<Ja)af3{<Jb)'^^ = -VabSc^ - {cTabjai 
{(^ajafsicrbcjj^ = -'>la\b{'^c]) a-f - i(^abc)aj, 

i<Ja)al3{<Jbcdf'' = -■^Va\bi(7cd])a " {<^abcd)<^, (9) 
{<^a)al3{crbcde)j^ = ■^Valbi'^cde]) ai + {(^abcde)aj, 
{(^a)al}{(^bcdef)'^'^ = " ^%[6(<^cde/] ) J " ^ea6cde/[4] (^'^^ ) J. 

{(7a)al3{(^b)'^^ = -VabSJ - ((Tab) J, 

M0''{<Taf'' = +Va\b{(rc]r - {(Tabor, 

{(rbcd)a0{(Taf'^ = --^Va\b{(Tcd\)a + {<Tabcd)^ , (10) 

{<Jbcde)p{<Ja)^^'' = -^??afb(q'cde])°^ + {( ^ abcde T^ , 

{<ybcdef)al3{<ya)^^ = " ^%[b((Tcde/] ) J + ^eabcde/[41 (o"'"' )c7- 

All of the sigma matrices are totally antisymmetric in their vector indices. The sigma matrices with 
one or five vector indices are symmetric with respect to spinor indices, while the matrix with three 
vector indices is antisymmetric in spinor indices. The symmetrization and antisymmetrization is 
denoted in the following way: 

= AaBp - AffBa- (11) 

The sigma matrices with five vector indices satisfy the identities: 



(a-[5])a/3 = 5je[5][5](<^ )a/3' 
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It is also useful to define matrices with more than five indices 

(o-[7])a/3 = -^e[7][3](CTt^l)a/3, 

i'^m).^ = (13) 

(o-[9])a/3 = e[9]a(0'-)a/3, 
\aB ^ ( „a\aB 



The trace identities are: 



{cJa)a0{<J^r^ = -16 4 ' 

= -16 45f<5^ , (14) 



/ c-S eft. c-i c-i . ll—h 

abode i 



In terms of ordinary matrices, our notation can be regarded as follows. We can introduce a bar- 
notation to distinguish the altitude of the indices on the bi-spinors in Eq(5). In otherwords we define 
the matrices 

CTa = (i7a)a/3 , = {cTa)"'^ , (15) 

and similarly for the third and fifth rank space-time tensors. We further use a Northwest-Southeast 
convention with respect to the contraction of hidden spinor indices. This is simplest to understand 
via the following examples 

<7-a,^6 = -(<7-a)a/3(<7-5)^^ ' ^b(^a = {(^bp'^{(^a)0a ■ (16) 

So when the contraction is up-to-down there is a plus sign when writing out the explicit indices, but 
when the contraction is down-to-up there is a minus sign. 

The sigma matrices with more vector indices follow from the multiplication table; 

(^apb = Vab + ab , 

^a<7b = -rjab + ab , 

(JafJbc — Valbf^c] + abc , 

CTbcCTa = -%[bCrc] -I- aabc , 
1 

(^a<^bcd = ^%[5<7cd] + O 'abcd , 

^bcd^a = 2[%[6'^cd] + O 'abcd , (17) 
_ 1 

CgO 'bcde — ~-^Va\b^cde] " O " abode ) 
4 



<^a ^bcde f_ 



1 

— yValbf^cde] + O 'abcde , 

1 ' ^1 [4] 

J^r}a\bfcder\ + ^eabcde/[4]0-^ 

1 1 



bcdelCTa— - ^Va\b<^cdel] - ^ ea6cde/[4] O" 



[4] 



CaC^fec = 


1la\b<^c] + f^abc , 


O'bc'^n = 






1 


O^ggfecd = 






1 


CTbcd^a = 






1 _ 


O'aP'bcde = 


-^Valbf^cde] + fabcde , 




1 


<^bcde^a = 


-^Va\b<^cde] — abode , 




1 1 1 


O'gCrbcde.f = 


--^,Va\b<^cdef] + ^ea6cde/[4]i7' 


<^bcdef^a = 


1 1 [4] 

-^Va\b(^cdef] - ^ea6cde/[4]0-' 



(18) 



[4] 



This establishes the enveloping algebra of the 10 D chiral matrices. It is easy to see that one can 
use these rules to reduce any product of such matrices down to a linear combination of irreducible 
elements. However the presence of the two copies of five forms down to one forms renders many 
of the calculatons intractable. In what follows we will introduce a program using the Mathematica 
language that can facilitate such computations. We will first give an example worksheet that shows 
how to use the program then present results that have exploited the routine. In the last section we 
present the program verbatim. We will call the program Sigma VectorlOD.m. 



II. INSTRUCTIONS FOR USING SIGMAVECTORIOD.M 



Here is an example run of the Mathematica package Sigma VectorlOD.m. First one must copy the 
program over to a Mathematica notebook and save the file as Sigma VectorlOD.m. Then open 
another notebook and bring in this file by typing, "<< Sigma Vector. m". Below is a transcript of 
a typical run. 

In[l]:= <<SiginaVectorlOD .m 
Out[l] = 

This is Sigma Vector program. It reduces the product of sigma matrices. Spinor indices of the 
chiral sigma matrices analyzed here are supressed. Sigma matrices with spinor indices down are 
denoted by S, as well as those matrices which have one spinor index down and other up. If 
both spinor indices are up, the letter SB is used. The program handles all contraction possi- 
bilities of spinor indices, except when two sigma matrices, both with even number of vector in- 
dices are multiplied in such a way that the lower index of the first sigma matrix is summed with 
the upper index of the second matrix. In this case the order of sigma matrices should be re- 
versed. Here is an example: {<j''-^)a * (cabc)"^ * (^e/)^ should be represented in the symbolic form 
as: Exp£indAll[ExpandAll[S[e, f] * *SB[a, b, c]] * *S[e, f]//.Rulel], where Rulel is a rule which 
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transforms sigma matrices with 6 or more vector indices in appropriate sigma matrices with 4 or 
less vector indices and contracts corresponding epsilon symbols. 

In [2]:= Here is the first example: 

In[3] := Expand[S[a, b] * *S[a, b, c] ] 

Out [3]= -2 S[c]+3 Ten S[c]-Ten^ S[c] 

In [4]:^ The output tells you something about the Kroneker 

delta symbols. In most cases you are not interested in 

that, so you type: 
In [5] := Ten = 10 
Out[5]= 10 

In[6}:= The previous example now becomes: 

In[7] := Expand[S[a, b] * *S[a, b, c] ] 

Out[7]= -72 S[c] 

In[8] := ExpandAll [ 

ExpandAll [S [e, f ] * *S [a] ] * *S [e, f ] ] 

Out[8]= 54 S[a] 

In [9]: = 

In [10] := Note that we have grouped the products of sigma matrices. 

This reduces the run time by a factor of three! When you 
have more than 10 vector indices, this is very important. 
In this simple example we could have used the following 

command : 

In[ll] := S[e, f] * *S[a] * *S[e, f] 

Out[ll]= 54 S[a] 

In [12]:= Here comes a nontrivial example: 

In[13] := ExpandAll [ 

ExpandAll [SB [e, f, g] * *S[a, b, c]] * * 
SB[e, f, g]] 

Out[13]= 258 SB[a, b, c]+^ Epsilon [a, b, c, e, 
f, g, $1[1] , $2[1] , $3 [1] , $4 [1] ] 
SB[e, f, g, $1 [1] , $2 [1] , $3 [1] , $4 [1] ] 
In [14]:= We have obtained a sigma matrix with 7 vector indices! 

The introduction of such matrices is important - the run 
time is much shorter. We now want to transform this 
7-index matrix into a product of epsilon symbol and a 
matrix with 3 indices by the use of "Rulel". Two epsilon 
symbols will be automatically contracted. 
In [15] := %//. Rulel 
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Out[15]= 258 SB[a,b,c]- 

(-5040 Id[a,$3[6]] Id[b,$2[6]] 

Id[c, $1 [6] ] + 5040 Id[a, $2[6]] 

Id[b, $3[6] ] Id[c, $1[6] ] + 

5040 Id[a, $3 [6] ] 

Id[b, $1 [6] ] Id[c, $2 [6] ] - 

5040 Id[a, $1[6]] Id[b, $3[6]] 

Id[c, $2 [6] ] - 5040 Id[a,$2[6]] 

Id[b, $1 [6] ] Id[c, $3 [6] ] + 

5040 Id[a, $1[6]] 

Id[b, $2 [6] ] Id[c, $3 [6] ] ) 

SB[$1 [6] , $2 [6] , $3 [6] ] 

In [16]:= In the previous input command line the % symbol stands 
for "the last output". Similarly, %% stands for the 
output before the last one. Better way is to use the 
output name directly, like in Out [ 17 ] / / . Rulel . 
The command //.Rulel means: apply Rulel as many times 
as it is applicable. If you want to apply it only once, 
the appropriate command is /.Rulel. 

What we have obtained in the last output is the correct 
answer in an unwanted form. To get the final answer we 
type: 

In[17] := ExpandAll [ % ] 

Out [17]= 48 SB [a, b, c] 

In [18] := When we put all this together, we see that the optimal 
command should have been: 

In[19] := ExpandAll [ 

ExpandAll [SB [e, f , g] * *S[a, b, c] ] * * 

SB[e, f, g]//. 

Rulel] 

Out [19]= 48 SB [a, b, c] 

In[20]:= It is always better to use ExpandAll then Expand. 

Note also that Rulel has been used only once and at the 

end of the calculation! If you use it more than once in 

the same command line you will get wrong answer! 

Sometimes, you want to use only a portion of the intermediate 

result and multiply it with another sigma matrix - to check 

the result, for example. 

In [21 ] := ExpandAll [SB [e, f , g] * *S[a, b, c] ] 
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Out [21] = -{Etala, g] Eta[b, f] Eta[c, e]) + 

Eta[a, f] Eta[b, g] Eta[c, e] + 

Eta[a,g] Eta[b, e] Eta[c, f]- 

Eta[a, e] Eta[b, g] Eta[c, f]- 

Eta[a, f] Eta[b, e] Eta[c, g] + 

Eta[a, e] Eta[b, f] Eta[c, g] — 

Eta[b, g] Eta[c, f] S[a, e] + 



Eta 


b. 


f ] 


Eta [c. 


g] 


S [a, e] + 


Eta 


b, 


g] 


Eta [c, 


e ' 


S[a, f]- 


Eta 


b, 


e ] 


Eta [c, 


g] 


S[a, f]- 


Eta 


b. 


f ] 


Eta[c, 


e 


S [a, g] + 


Eta 


b. 


e] 


Eta [c. 


f ] 


S[a, g] + 


Eta 


a. 


g] 


Eta [c. 


f ] 


S [b, e] - 


Eta 


a. 


f ] 


Eta [c. 


g] 


S [b, e] - 


Eta 


a. 


g] 


Eta [c. 


e; 


S[b, f] + 


Eta 


a. 


e] 


Eta [c. 


g] 


S [b, f ] + 


Eta 


a. 


f ] 


Eta[c, 


e 


S[b, g]- 


Eta 


a, 


e ] 


Eta[c, 


f] 


S [b, g] - 


Eta 


a, 


g] 


Eta [b, 


f] 


S [c, e] + 


Eta 


a, 


f ] 


Eta[b, 


g] 


S [c, e] + 


Eta 


a. 


g] 


Eta[b, 


e; 


S[c, f]- 


Eta 


a. 


e] 


Eta[b, 


g] 


S[c, f]- 


Eta 


a. 


f ] 


Eta[b, 


e; 


S[c, g] + 


Eta 


a. 


e] 


Eta[b, 


f ] 


S[c, g]- 


Eta 


c, 


g] 


S [a, b, 


e. 


f] + 


Eta 


c, 


f ] 


S [a, b, 


e. 


g]- 


Eta 


c. 


e ] 


S [a, b. 


f , 


g] + 


Eta 


[b. 


g] 


S [a, c. 


e. 


f]- 


Eta 


b. 


f ] 


S [a, c. 


e. 


g] + 


Eta 


b. 


e] 


S [a, c. 


f , 


g]- 


Eta 


a. 


g] 


S [b, c. 


e. 


f] + 


Eta 


a, 


f ] 


S [b, c, 


e. 


g]- 


Eta 


a. 


e ] 


S [b, c, 


f , 


g]- 


S [a 


b. 


c. 


e, f , g] 







In[22]:= Eta symbol is the same as Id - depending whether both 
vector indices are down or up or not. Suppose now that 
you want to extract the last term S [a, b, c, e, f , g] and 
to multiply it with SB[e,f,g] . You first click the mouse 
on the last output and press "Command-x", i.e. the 
Command key followed by x. This cuts the output. Put it 
back by pressing "Command-v". This will paste it back. 
Move the cursor beneath the last cell and press again 



"Command-v" . 
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Out [22] = -{Etala, g] Eta[b, f] Eta[c, e]) + 



Eta 


a. 


f ] 


Eta[b, 


g] 


Eta[c, 


e] + 


Eta 


a. 


g] 


Eta [b, 


e ' 


Eta 


c, 


f]- 


Eta 




e ] 


Eta [b, 


g] 


Eta 


c, 


f ] - 


Eta 


a. 


f ] 


Eta [b, 


e 


Eta 


c, 


g] + 


Eta 


a. 


e] 


Eta[b, 


f ] 


Eta[c, 


g]- 


Eta 


[b. 


g] 


Eta [c. 


f ] 


S [a. 


e 


+ 


Eta 


[b. 


f ] 


Eta [c. 


g] 


S [a. 


e 


+ 


Eta 


[b. 


g] 


Eta [c. 




S [a. 


f 


— 


Eta 


b. 


e] 


Eta [c, 


g] 


S [a. 


f 


— 


Eta 


b. 


f ] 


Eta [c, 


e 


S [a. 


g] 


+ 


Eta 


b, 


e ] 


Eta [c, 


f ] 


S [ a. 


g] 


+ 


Eta 


a. 


g] 


Eta [c, 


f ] 


S [b. 


e 


— 


Eta 


a. 


f ] 


Eta [c. 


g] 


S [b. 


e 


— 


Eta 


a. 


g] 


Eta [c. 


e; 


S [b. 


f 


+ 


Eta 


a. 


e] 


Eta [c. 


g] 


S [b. 


f 


+ 


Eta 


a. 


f ] 


Eta [c. 


e; 


S [b. 


g] 


— 


Eta 


a. 


e ] 


Eta [c, 


f ] 


S [b. 


g] 


— 


Eta 


a. 


g] 


Eta [b, 


f ] 


S [ c. 


e 


+ 


Eta 


a. 


f ] 


Eta [b. 


g] 


S [c. 


e 


+ 


Eta 


a. 


g] 


Eta[b, 


e; 


S [c. 


f 


— 


Eta 


a. 


e] 


Eta [b. 


g] 


S [c. 


f 


— 


Eta 


a. 


f ] 


Eta[b, 


e; 


S [c. 


g] 


+ 


Eta 


a. 


e] 


Eta [b. 


f ] 


S [c. 


g] 


— 


Eta 


c. 


g] 


S[a, b. 


e. 


f] + 






Eta 


c. 


f ] 


S [a, b, 




g]- 






Eta 


c, 


e ] 


S [a, b, 


f, 


g] + 






Eta 


b. 


g] 


S [a, c, 


s, 


f]- 






Eta 


[b. 


f ] 


S [a, c. 


e. 


g] + 






Eta 


[b. 


e] 


S [a, c. 


f , 


g]- 






Eta 


a. 


g] 


S[b, c. 


e. 


f] + 






Eta 


a. 


f ] 


S[b, c. 


e. 


g]- 






Eta 


a. 


e] 


S [b, c, 


f , 


g]- 






S [a 


b 


c. 


e, f, g] 











In[23]:= Now you want this output to be turned into an input - 

bold face. Click again on it and press "Command 9". This 
will turn it into a bold-faced form: 
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In[24] := —{Etala, g] Eta[b, f] Eta[c, e]) + 



Eta[a, f] 


Eta[b, 


g] 


Eta[c, e]-f 


Eta [a, g] 


Eta[b, 


e] 


Eta[c, f ] — 


Eta [a, e] 


Eta[b, 


g] 


Eta[c, f ] — 


Eta[a, f] 


Eta[b, 


e] 


Eta[c, g] + 


Eta[a, e] 


Eta[b, 


f ] 


Eta[c, g] — 


Eta[b, g] 


Eta[c, 


f ] 


S[a, 


e] + 


Eta[b, f] 


Eta[c, 


g] 


S[a, 


e] + 


Eta[b, g] 


Eta[c, 


e] 


S[a, 


f ] — 


Eta[b, e] 


Eta[c, 


g] 


S [a. 


f ] — 


Eta[b, f] 


Eta[c, 


e] 


S [a. 


g] + 


Eta[b, e] 


Eta[c, 


f ] 


S [a. 


g] + 


Eta [a, g] 


Eta[c, 


f ] 


S [b. 


e] — 


Eta[a, f] 


Eta[c, 


g] 


S[b, 


e] — 


Eta [a, g] 


Eta[c, 


e] 


S[b, 


f ] + 


Eta[a, e] 


Eta[c, 


g] 


S[b, 


f ] + 


Eta[a, f] 


Eta[c, 


e] 


S[b, 


g] — 


Eta [a, e] 


Eta[c, 


f ] 


S [b. 


g] — 


Eta [a, g] 


Eta[b, 


f ] 


S [c. 


e] + 


Eta [a, f] 


Eta[b, 


g] 


S[o, 


e] + 


Eta [a, g] 


Eta[b, 


e] 


S[c, 


f ] — 


Eta[a, e] 


Eta[b, 


g] 


S[c, 


f ] — 


Eta[a, f] 


Eta[b, 


e] 


S[c, 


g] + 


Eta[a, e] 


Eta[b, 


f ] 


S[c, 


g]- 


Eta[c, g] 


S[a, b, 


e, 


f] + 




Eta[c, f] 


S [a, b. 


e. 


g]- 




Eta [o, e] 


S [a, b. 


f , 


g] + 




Eta[b, g] 


S [a, c, 


e, 


f]- 




Eta[b, f] 


S [a, c, 


e. 


g] + 




Eta[b, e] 


S[a, c, 


f, g]- 




Eta [a, g] 


S[b, c, 


e, 


f] + 




Eta[a, f] 


S[b, c, 


e, g] — 




Eta[a, e] 


S[b, c, 


f, g]- 




S [a, b, c. 


e, f , g] 









In [25]:= You can now darken the part which you want to throw out: 
Click the mouse at the beginning of the last input and 
drag it until the last term. Cut the darkened piece 
by "Command x" . What is left is S [a, b, c, e, f , g] . 

In [26]:= Below are two examples that correct Eqs(40) and (43) in the 
appendix of [6] : 

In[27] := ExpandAll [ SB [ f , a, b] * *S[f , c, d] ] 
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Out[27]= -8 Eta[a, d] Eta[b, c] + 

8 Eta[a, c] Eta[b, d]-7 Eta[b, d] S[a, c] + 
7 Eta[b, c] S[a, d]+7 Eta[a, d] S[b,c]- 
7 Eta[a, c] S[b, d]-6S[a,b, c, d] 

In [28] := ExpandAll [SB [a, b, c, d, e] * *S[a, f , g] ] 

Out[28]= 6 Eta[d, g] Eta[e, f] S[b, c]- 



6 


Eta[d, f ] 


Eta[e, 


g] 


S[b, c]- 


6 


Eta [c, g] 


Eta[e, 


f ] 


S [b, d] + 


6 


Eta[c, f] 


Eta [e. 


g] 


S [b, d] + 


6 


Eta [c, g] 


Eta[d, f] 


S [b, e] - 


6 


Eta[c, f] 


Eta [d, g] 


S [b, e] + 


6 


Eta [b, g] 


Eta [e. 


f ] 


S [c, d] - 


6 


Eta[b, f ] 


Eta[e, 


g] 


S[c, d]- 


6 


Eta[b, g] 


Eta[d, f] 


S [c, e] + 


6 


Eta[b, f ] 


Eta [d, g; 


S [c, e] + 


6 


Eta[b, g] 


Eta[c, 


f ] 


S[d, e]- 


6 


Eta [b, f ] 


Eta [c. 


g] 


S [d, e] + 


5 


Eta [e, g] 


S [b, c. 


d. 


f]- 


5 


Eta[e, f] 


S [b, c. 


d. 


g]- 


5 


Eta[d, g] 


S [b, c. 


e. 


f] + 


5 


Eta[d, f ] 


S [b, c. 


e. 


g] + 


5 


Eta [c, g] 


S [b, d. 


e. 


f]- 


5 


Eta[c, f] 


S [b, d. 


e. 


g]- 


5 


Eta[b, g] 


S[c, d. 


e. 


f] + 


5 


Eta [b, f ] 


S [c, d. 


e. 


g] + 


4 


S [b, c, d. 


e, f, g] 







III. IDENTITIES AND CALCULATED RESULTS 

Before using the program to find new identities we use the definitions from the first section to 
construct the following Fiertz identities. Note that our program will not derive these identities since 
they are not products of the elements via matrix multiplication. 



{<y-)ic.(}{cTa)y)S = 0, (19) 

{'^-)ail3{<7a)j)S = -(CT-)/37(CTa)a5, (20) 

{a^)'^(P{a^y)' =-{a^)P^aar\ (21) 

(^''l)a/3(cT[5])7i = 0, (22) 

{cT^)a0{(Tabcr' = S ■ 31 Sjjf,^ , (23) 

(cr—)a/3 ((7060)75 = -2 • 3! {(T-)o:['y{o'a)s]0, (24) 

(^£6c)a/3(^^)75 = _2 . 3! (a^)«b(^^)*l/3, (25) 

(^''') J(^[4])/ = 4!{-2 5ZS'^ + 12 SjSi - 2 {a^U{a^r% (26) 
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(^['l)aMq5])^' = 5!{-16 - 2 {a^Uia^r'}, (27) 

(a^) J(cTa6)/ = 2 {-5l5l - 4 5^^^ - 2 (a^)„;3K)^n, (28) 

= &5U%, (29) 

= - 8 ('^-)«/3('^a)^', (30) 

[2Sjj'0 + {a'-\^f,\{a,y']{a^Us = 0, (31) 

{{<^\a)Hj{cr'-^r' + l{<r'-^)H-y{<riar'K'^bc^h)5 = 48 {a^^U ^f^^' (^2) 

('7'a6c)(a(/3(<7'-)7)5) = "2 (o-[o)c<5(ct6])/3t,, (33) 

(0--)(a/3(CTa5)^') = - (t7'&)(a/35^) , (34) 

(a^)("''(a^)/^ = (<76)("''^J\ (35) 

(0'-)(a/3('^a6)-y')(CTc)6e = -2(a-[a)(a/3(cr6])-y)e. (36) 

There is also an interesting identity for which it is convenient to first introduce two arbitrary three- 
forms in order to write the identity in a simple manner, 

(crl^l— )(£,^|(cr[^ lfe;)|^5)A[3]i3[3/] = (a--)(a/3| {a-)\^s)^acdBb— 

- 6{a-)(ai3\{a^^-)\^s){AabcBdef^ + BgbcAdef) , (37) 

valid for the arbitrary A[3] = Aabc and S[3'] ^ Bdef- 

Finally, the following identities for manipulating the cr-matrices are consequence of our definitions 
and are a direct test of the reduction program, Sigma VectorlOD.m. 

{<JbcdU{<J^f^ = 72 6fSZ + 56(a/)J, (38) 

{a^)J'{ab_cd)0y = -7 • 8 • 9 {a^)^^, (39) 

^a^)^0{a^)p, = 8 {a^\", (40) 

(a^)a/3(ae)^^ = -8 (a^)J, (41) 

(aa)aMcT^)/ = -9 {a'U, (42) 

i'Jar^ia^); = 9 iaT', (43) 

[cT^'X^iaah, = -7 (a^)„„ (44) 

(^a5cd)^a(^^)/37 = _7 ((j6£d)«7^ (45) 

(aa6)/(a^^)/ = - S^S^SZ - S^Ja/)J + {aJ^)J, (46) 

{aa_b,U(-^lr = sf^S^S^SZ + l6^Sl{a/v " l^^i-bj-^ " (- (47) 

(CT^)i'(o-afcc)/J7 - -72(a,)„,,, (48) 

{<Jabcr^ia^); = -72 (a,)"', (49) 

{a^r^iaabjp^ = -6-7ia^r\ (50) 

ic7^)Maab\>' = 6-7{a^)aj, (51) 

(^^)/(^a6-V = -7-8(a^)J, (52) 

(^-)/3"('^a.-)7^ = -7-8ia^)-, (53) 

(a^)"'^(aa6c)0^ = -6-7-8 (ai^)^«, (54) 

(a^)«'^(ahede)/ = 6 • 7 • 8 ■ 9 (a^)"T, (55) 

12 



a^r^{aa)f3^ = -6 

<T-cr'^{aabi)a^ = 7-8 {acd\^ - 8 • 9vciS^, 

a^)"'^(aa6/)^^ = -5 ■ 6 ((7^—)^" - ^4'(^'^'')7"' 

ae-)"''(<Tad)/(aefc)„^' = 3-19 {a^y if c^d^e, 



a6 



^[3])"''M/(Ma = ^^a6cd[3][3]('^'''r' [3] 7^ 6, C, ^ , 

0'-)(a/3(0-a6cd)^') = 2 (^[i) ("/3 (^^bc] )-/) , 

<''[a)(a/3(0'6c])-y'j = ~ fag (o"-)*^^ (o"abc)7)e > 

'7-)/((Tc)/37(<^a6)/ = 54(ae)a5, 

^-)/(^-)/(^a.)/ = -26 
aae)g"(a^^)/ = -8 (a^e)^" - 9 jfc^;^, 

crMcT^f^ = -WSZ, 

= -10-9 52, 

(Tabc)M(T"-^f^ = 10-9-8 53;, 

cTa6cd)/(a^)/ = 10-9-8-7 52, 

^abcde)ag(a^)^^ = 10 • 9 • 8 ■ 7 ■ 6 , 
q--)a/3 (q-gbc)^^ (o"/)7'5 = 4 (cro6c)c<5, 
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{(T^^-^)£{(Tahj{(Tl^lJs' = 54 {aa)a5, 

{^^^^^).^M,%<^tjJs' = -26 (^06)/, 

((7^1-2 )/((Ta6c)/37(f^/^/^)/ = 6 {(Tabc)a6, 

{c7U.)£[aa^Y'<[af_jJ^' = 6 ia^f\ 

(a--i-2)J'((Tabcde)/37(0'/^/^)/ = -10 (cra6cde)a5, 

= -10 {agbedef^ 
= -48 

('7-i-2-3)a/3(CTafcc)^'*'(cr/^/^/^)75 = 48 {(Tabc)a6, 

{aLLLrP(a^)^,^{a^j_^f_f' = 48 
(aZi/2/3)„^(a^)/(a/^/^/^)'^^ = -48 (aabed)a , 

^aLLLL)0{a^r{a^^^^^^^J,'' = 24-42 
(^/././3/.)/(a^)^«(a^^^^^^^J^^ = -24-14 (a^)/, 
{aLLLL)0{^j^)^^(af_j_j_j_^)^'< = -24-14 
(^/./.Z3£4)/(a^)«7(^^^^^^^^)^^ = -24-14(a^)''^ 

{aLiLLL)£{aabcde)l3^{cTf_J_^f_J_^)^ = 240 (CTa6cde)a5, 



IV. THE PROGRAM 

(*BeginPackage["SigniaVector"']*) 

Print ["This is Sigma Vector program. It reduces the product of sigma matrices. 
Spinor indices of the chiral sigma matrices analyzed here are supressed. 
Sigma matrices with spinor indices down are denoted by S, as well as those 
matrices which have one spinor index down and other up. If both spinor indices 
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are up, the letter SB is used. 

The program handles all contraction possibilities of spinor indices, except when 
two sigma matrices, both with even number of vector indices are multiplied in such 
a way that the lower index of the first sigma matrix is summed with the upper index 
of the second matrix. In this case the order of sigma matrices should be reversed. 
Here is an example: 

/})_{\a}"{\6} * (\,s.{«. h, c}y {\a\g} * (^e, 
should be represented in the symbolic form as: 
ExpandAll[ExpandAll[S[e, /] * *SB[a, b, c]] * *S[e, f]//.Rulel], 

where Rulel is a rule which transforms sigma matrices with 6 or more vector indices 
in appropriate sigma matrices with 4 or less vector indices and contracts corresponding 
epsilon symbols."] 

(* These next rules assumed repeated indices for implied sum *) 

Set Attributes [Id,Orderless] 

Ten=10 

Set Attributes [Eta, Orderless] 

Eta/: Eta[a_,b_]"2:= Ten (* 10 for our purposes *) 
Eta/: Eta[a_,a_]:= Ten 
Eta/: Eta[a_, b.]*S[d...,a_, c.._] := S[d,b, c] 
Eta/: S[d...,a_, c.__]*Eta[a., b_] := S[d,b, c] 
Eta/: Eta[a_, b.]*SB[d___,a_, c_..] := SB[d,b, c] 
Eta/: SB[d__.,a_, c_..]*Eta[a_, b.] := SB[d,b, c] 

Eta/: Eta[a_, b_]*Epsilon[d ,a_, c ] := Epsilon[d,b, c] 

Eta/: Epsilon[d ,a_, c ]*Eta[a_, b_] := Epsilon[d,b, c] 

Eta/: Eta[a_, b.]*Eta[a_, c_] := Eta[b, c] 

Eta/: Eta[a_, c_]*Eta[a_, b_] := Eta[b, c] 

Eta/: NumberQ[Eta[x__]] := True 

(* These rules define Ten as the dimension so that its 

dependence in the 

calculations can be followed. *) 

Unprotect [NumberQ] 

NumberQ [Ten ' n Jnteger] : =True 

NumberQ[Ten] :=True 
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NumberQ[l/four!]:=True (* used to denote 1/4! *) 
NumberQ[l/three!]:=True (* used to denote 1/3! *) 
Protect [Number Q] 

(* Set the tag parameter for implied sums here *) 
epi = 1 

resct:=:(cpi=l) 

Par [exp__] : = (oxp) 

(* This comes from Epsilon.m *) 

(* Rules for Epsilon when Euclidean space metric is present *) 
(* These Rules are true for any dimension. Just change Ten to 
the appropriate dimension. See also Rule2. *) 
Epsilon/ : NumberQ[Epsilon[x__]]:=True 

Epsilon/: Epsilon[b..]*Epsilon[b_.] := Ten!/; Signature[{b}]=!=0 
Epsilon/: Epsilon [b..]~ 2 := Ten! 
Unprotect[NonCommutativeMultiply] 

(x.?NumbcrQ a_[i..])**(y_?NumberQ b_[j__]):=(x y) a[i]**b[j] 
(a_ +b_)**c_ := a**c +b**c 

a. **(b_ + c_) := a**b + a**c 
a_ **(x_?NumberQ b_):= x a**b 
(x.?NumberQ a_)**b. := x a**b 

b. ?NumberQ ** a_ := b a 
a. ** b.?NumberQ := b a 
((anything__)*c_)**d_ := (anything) (c**d) 
Protect [NonCommutativeMultiply] 

(* Definition of an antisymmetrizer. *) 
SetAttributes[AntiSymmetrize,HoldAll] 
AntiSymmetrizc/: AntiSymmctrize[ k__,exp_] := 
Sum [Signature [Permutations [k] [ [$i]] ] *Signaturc [k] * 
(exp /. Table[k[[$j]] - > Permutations[k][[$i]][[$j]], 
{$j, Length[k]}]), 
{$i, 1, Length[k]!}] 
(* These rules define S,SB and Epsilon as purely 
anti-symmetric tensors 
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for any rank. *) 

S[x__] := Signature[{x}] S@@Sort[{x}] /; Not[OrderedQ[{x}]] 

SB[x..] := Signature[{x}] SB@@Sort[{x}] /; Not[OrderedQ[{x}]] 

Epsilon[x__] := Signature [{x}] Epsilon@@Sort[{x}] /; 

Not[OrdercdQ[{x}]] 

S[x__]:=0 /; Signaturc[{x}]==0 

SB[x..]:=0 /; Signaturc[{x}]==0 

Epsilon[x__]:=0 /; Signaturc[{x}]==0 

(* Epsilon[l,2,3,4,5,6,7,8,9,10] = 1 *) 

(* This rule replaces 6,7,8,9, and 10 forms with their duals *) 
Rulel = { 
S[a_,b_,c_,d_,e_,f_] :> 

( Epsilon[a,b,c,d,e,f,$l[epi],$2[epi],$3[epi],$4[epi]]* 
S[$l[epi],$2[epi],$3[epi],$4[epi++]])/24, 

S[d[l],d[2],d[3],d[4],d[5],d[6]]:> 

(Epsilon[d[l],d[2],d[3],d[4],d[5],d[6],$l[epi],$2[epi],$3[epi],$4[epi]]* 

S[$l[epi],$2[epi],$3[epi],$4[epi++]])/24, 

SB[a_,b_,c_,d_,e_,f_,g_] :> 

(- Epsilon[a,b,c,d,e,f,g,$l[epi],$2[epi],$3[epi]]* 

SB[$l[epi],$2[epi],$3[epi++]])/6, 

SB[d[l],d[2],d[3],d[4],d[5],d[6],d[7]]:> 

(-Epsilon[d[l],d[2],d[3],d[4],d[5],d[6],d[7],l[epi],$2[epi],$3[epi]]* 

SB[$l[epi],$2[epi],$3[epi++]])/6, 

S[a_,b_,c_,d_,e_,f_,g_] :> 

(- Epsilon[a,b,c,d,e,f,g,$l[epi],$2[epi],$3[epi]]* 

S[$l[epi],$2[epi],$3[epi++]])/6, 
S[d[l],d[2],d[3],d[4],d[5],d[6],d[7]]:> 

(-Epsilon[d[l],d[2],d[3],d[4],d[5],d[6],d[7],$l[epi],$2[epi],$3[epi]]* 
S[$l[epi],$2[epi],$3[epi++]])/6, 

S[a_,b_,c_,d_,e_,f_,g_,h_] :> 

(-Epsilon[a,b,c,d,e,f,g,h,$l[epi],$2[epi]]* 

S[$l[epi],$2[epi++]])/2, 
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S[d[l] ,d[2] ,d[3] ,d[4] ,d[5] ,d[6] ,d[7] ,d[8]] :> 

(-Epsilon[d[l],d[2],d[3],d[4],d[5],d[6],d[7],d[8],$l[epi],$2[epi]]* 

S[$l[epi],$2[epi++]])/2, 

SB [a_,b_,c_,d_,e_,f_,g_,h_,i_] : > 

(Epsilon[a,b,c,d,e,f,g,h,i,$l[epi]]* 
SB[$l[epi++]]), 

SB[d[l],d[2],d[3],d[4],d[5],d[6],d[7],d[8],d[9]] :> 

(Epsilon[d[l],d[2],d[3],d[4],d[5],d[6],d[7],d[8],d[9],$l[epi]]* 

SB[$l[epi++]]), 

S[a_,b_,c_,d_,e_,f_,g_,h_,i_] :> 

(Epsiloii[a,b,c,d,e,f,g,h,i,$l[epi]]* 

S[$l[epi++]]) 

} 

(* This rule pulls out the coefficients of chiral matrices. *) 
Set Attributes [Eta2 ,Orderless] 
Eta2/: NumberQ[Eta2[x_]] := True 
Rulc2 = { 

S[a_] :> S[d[l]] * Eta2[d[l],a], 

S[a_,b.] :> S[d[l],d[2]] * Eta2[d[l],a] Eta2[d[2],b], 

S[a_,b.,c.] :> S[d[l],d[2],d[3]]*Eta2[d[l],a]* Eta2[d[2],b] Eta2[d[3],c], 

S[a_,b.,c.,dl_] :> S[d[l],d[2],d[3],d[4]] * 

Eta2[d[l],a]* Eta2[d[2],b] Eta2[d[3],c] Eta2[d[4],dl], 

S[a_,b.,c.,dl_,e.] :>S[d[l],d[2],d[3],d[4],d[5]] * 

Eta2[d[l],a]* Eta2[d[2],b] Eta2[d[3],c] Eta2[d[4],dl] Eta2[d[5],e], 

S[a_,b.,c.,dl-,e-,f-] :>S[d[l],d[2],d[3],d[4],d[5],d[6]] * 

Eta2[d[l],a]* Eta2[d[2],b] Eta2[d[3],c] Eta2[d[4],dl] Eta2[d[5],e] * 

Eta2[f,d[6]], 

S[a_,b.,c.,dl_,c.,L,g_] :>S[d[l],d[2],d[3],d[4],d[5],d[6],d[7]] * 
Eta2[d[l],a]* Eta2[d[2],b] Eta2[d[3],c] Eta2[d[4],dl] Eta2[d[5],e] * 
Eta2[f,d[6]] Eta2[g,d[7]], 

S[a_,b.,c.,dl_,e.,L,g_,h.] :>S[d[l],d[2],d[3],d[4],d[5],d[6],d[7],d[8]] * 
Eta2[d[l],a]* Eta2[d[2],b] Eta2[d[3],c] Eta2[d[4],dl] Eta2[d[5],e] * 
Eta2[f,d[6]] Eta2[g,d[7]] Eta2[h,d[8]], 
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S[a_,b.,c.,dl_,e.,L,g_,h.,L] :>S[d[l],d[2],d[3],d[4],d[5],d[6],d[7],d[8],d[9]] * 
Eta2[d[l],a]* Eta2[d[2],b] Eta2[d[3],c] Eta2[d[4],dl] Eta2[d[5],e] * 
Eta2[f,d[6]] Eta2[g,d[7]] Eta2[h,d[8]] Eta2[i,d[9]], 

S[a_,b.,c.,dl_,e.,f.,g_,h.,i.,j.] :>S[d[l],d[2],d[3],d[4],d[5],d[6],d[7],d[8],d[9],d[10]] * 

Eta2[d[l],a]* Eta2[d[2],b] Eta2[d[3],c] Eta2[d[4],dl] Eta2[d[5],e] * 
Eta2[f,d[6]] Eta2[g,d[7]] Eta2[h,d[8]] Eta2[i,d[9]] Eta2[j,d[10]], 
SB[a.] :> SB[d[l]] * Eta2[d[l],a], 

SB[a.,b.,c.] :> SB[d[l],d[2],d[3]]*Eta2[d[l],a]* Eta2[d[2],b] Eta2[d[3],c], 
SB[a.,b.,c.,dl_,e_] :>SB[d[l],d[2],d[3],d[4],d[5]] * 
Eta2[d[l],a]* Eta2[d[2],b] Eta2[d[3],c] Eta2[d[4],dl] Eta2[d[5],e], 
SB[a.,b.,c.,dl_,e_,f_,g.] :>SB[d[l],d[2],d[3],d[4],d[5],d[6],d[7]] * 
Eta2[d[l],a]* Eta2[d[2],b] Eta2[d[3],c] Eta2[d[4],dl] Eta2[d[5],e] * 
Eta2[f,d[6]] Eta2[g,d[7]], 

SB[a.,b.,c.,dl_,e_,f_,g.,h_,i_] :>SB[d[l],d[2],d[3],d[4],d[5],d[6],d[7],d[8],d[9]] * 
Eta2[d[l],a]* Eta2[d[2],b] Eta2[d[3],c] Eta2[d[4],dl] Eta2[d[5],e] * 
Eta2[f,d[6]] Eta2[g,d[7]] Eta2[h,d[8]] Eta2[i,d[9]] 
} 

(* Use Recover to factor out the S's and SB's *) 

Recover[exp_] := Collect [exp/.Rule2,{ S[d[l]], S[d[l],d[2]],S[d[l],d[2],d[3]], 

S[d[l] ,d[2] ,d[3] ,d[4]] , S[d[l] ,d[2] ,d[3] ,d[4] ,d[5]] , 

S[d[l],d[2],d[3],d[4],d[5],d[6]], 

S[d[l],d[2],d[3],d[4],d[5],d[6],d[7]], 

S[d[l] ,d[2] ,d[3] ,d[4] ,d[5] ,d[6] ,d[7] ,d[8]] , 

S[d[l],d[2],d[3],d[4],d[5],d[6],d[7],d[8],d[9]], 

S[d[l],d[2],d[3],d[4],d[5],d[6],d[7],d[8],d[9],d[10]], 

SB[d[l]], SB[d[l],d[2],d[3]], SB[d[l],d[2],d[3],d[4],d[5]], 

SB[d[l],d[2],d[3],d[4],d[5],d[6],d[7]], 

SB[d[l],d[2],d[3],d[4],d[5],d[6],d[7],d[8],d[9]]}]/.Rule3 

Rule3={Eta2[f_]:> Eta[f]} 

(* This rule is time consuming. *) 
Epsilon/: Epsilon[b__]*Epsilon[c__] := 
-Par[ Signature[ Join[Complement[{b}, 
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Intersection[{b},{c}]], 
Intersection[{b},{c}]]]* 

Signature[ Join[Complenient[{c},Intersection[{b},{c}]], 
Intersection[{b},{c}]]]*ReleaseHold[ ( 

AntiSymmetrize[Complenient[{c}, Intersection[{b}, {c}]], 
Product [Eta[Coinplcmcnt[{b}, Intcrscctioii[{b}, {c}]][[$i]], 
Complement [{c}, Intcrsection[{b}, {c}]][[$i]]], 
{$i, Length [Complement [{c}, Intersection[{b}, 
{c}]]]}]])]*(Length[ 

Intersection[{b}, {c}]]!)] /; {b}=!={c} && Length[{b}]==Length[{c}] 
S/: S[L]**SB[j_] := - Eta[i, j] - S[i, j] 
S/: S[i.]**S[j.,k.] := -S[k] Eta[i, j] + S[j] Eta[i, k] - 
S[i, j, k] 

S/: S[j.,k.]**S[i_] := -S[k] Eta[i, j] + S[j] Eta[i, k] + 

S[i, j, k] 

S/: S[j.,k.]**SB[L] := Eta[i,j] SB[k] - Eta[i,k] SB[j] - 
SB[i,j,k] 

S/: S[L] ** SB[j_, k., L] := 
-S[k, 1] Eta[i, j] - S[l,j] Eta[i, k] - S[j, k] Eta[i, 1] - 
S[i, j, k, 1] 

S/: Sp.,k.,L]**SB[i.] := -Eta[i,j] S[k,l] - Eta[i,k] S[l,j] - 

Eta[i,l] S[j,k] + S[i,j,k,l] 

S/: S[i.]**S[j.,k.,L,m.] := Eta[i,j] S[k,l,m] - 

Eta[i,k] S[l,m,j] + Eta[i,l] Sp,k,m] - 

Eta[i,m] S[j,k,l] + S[i,j,k,l,m] 

S/: S[j.,k.,L,m.]**S[L] := -Eta[i,j] S[k,l,m] + 

Eta[i,k] S[l,m,j] - Eta[i,l] S[j,k,m] + 

Eta[i,m] S[j,k,l] + S[i,j,k,l,m] 

S/: S[j.,k.,L,m.]**SB[L] := -Eta[i,j] SB[k,l,m] + 

Eta[i,k] SB[l,m,j] - Eta[i,l] SB[j,k,m] + 

Eta[i,m] SBp,k,l] + SB[i,j,k,l,m] 

S/: S[i.]**SBp_,k_,L,m_,n.] := -Eta[i,j] S[k,l,m,n] - 

Eta[i,k] S[l,m,n,j] - Eta[i,l] S[m,n,j,k] - 
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Eta[i,m] S[n,j,k,l] - Eta[i,ii] S[j,k,l,m] - S[i,j,k,l,m,n] 
S/: Sp.,k.,L,ni.,n_]**SB[L] := -Eta[i,j] S[k,l,ni,n] - 
Eta[i,k] S[l,ni,n,j] - Eta[i,l] S[m,n,j,k] - 
Eta[i,m] S[ii,j,k,l] - Eta[i,ii] S[j,k,l,m] + S[i,j,k,l,m,n] 

S/: S[L]**S[j.,k.,L,m.,n_,p.] := - S[i„i,k,l,m,n,p] + 
(Epsilon[j,k,l,ni,n,p,$l[cpi], $2[epi], $3[cpi], $4[epi]]* 
S[i,$l[epi], $2[cpi], $3[cpi], $4[epi++]])/24 
S/: S[j.,k.,L,m.,n_,p.]**SB[L] := SB[i,j,k,l,m,n,p] + 
(Epsilon[j,k,l,m,n,p,$l[epi], $2[epi], $3[epi], $4[epi]]* 
SB[i,$l[epi], $2[epi], $3[epi], $4[epi++]])/24 
S/: Sp.,k.,L,ni.,n_,p.]**S[L] := S[i,j,k,l,m,n,p] + 
(Epsiloii[j,k,l,m,n,p,$l[epi], $2[epi], $3[epi], $4[epi]]* 
S[i,$l[epi], $2[epi], $3[epi], $4[epi++]])/24 
S/: S[i.]**SB[j_,k_,l_,m_,n.,p_,q.] := S[i,j,k,l,ni,n,p,q] + 
(Epsilon[j,k,l,m,n,p,q,$l[epi], $2[cpi], $3[cpi]]* 
S[i,$l[epi], $2[epi], $3[epi++]])/6 

S/: S[j.,k_,L,m_,n_,p_,q_]**SB[L] := S[i,j,k,l,m,n,p,q] - 
(Epsilon[j,k,l,m,n,p,q,$l[epi], $2[cpi], $3[epi]]* 
S[i,$l[epi], $2[epi], $3[epi++]])/6 
S/: S[i.]**S[j.,k.,L,ni.,n_,o.,p.,q.] := 
S[i,j,k,l,m,n,o,p,q]+ 

(Epsiloii[j,k,l,m,ii,o,p,q,$l[epi], $2[epi]]* 

S[i,$l[epi], $2[epi++]])/2 

S/: S[j.,k.,L,ni.,n_,o.,p.,q_]**S[i.] := 

S[i,j,k,l,m,ii,o,p,q]- 

(Epsilon[j,k,l,m,n,o,p,q,$l[epi], $2[epi]]* 

S[i,$l[epi], $2[epi++]])/2 

S/: S[j.,k.,L,m.,n_,o.,p.,q_]**SB[L] := 

- SB[i,j,k,l,m,n,o,p,q]+ 

(Epsilon[j,k,l,m,n,o,p,q,$l[epi], $2[epi]]* 

SB[i,$l[epi], $2[epi++]])/2 

S/: S[i.]**SB[j_,k_,L,m_,n.,o_,p_,q.,r.] := 

Epsilon[i,j,k,l,m,n,o,p,q,r] - 
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Epsilon[j,k,l,m,n,o,p,q,r,$l[epi]]*S[i,$l[epi++]] 
S/: Sp.,k.,L,ni.,n_,o.,p.,q_,r_]**SB[L] := 
Epsiloii[i,j,k,l,m,ii,o,p,q,r] + 
Epsiloii[j,k,l,m,ii,o,p,q,r,$l[epi]]*S[i,$l[epi++]] 
SB/: SB[L] ** S[j.] := -Eta[i, j] + S[i, j] 
SB/: SB[L]**S[j_,k.] := Eta[i,j] SB[k] - Eta[i,k] SB[j] + 
SB[i,j,k] 

SB/: SB[L]**S[j_,k.,L] := Eta[i,j] S[k,l] + 
Eta[i,k] S[l,j] + Eta[i,l] S[j,k] - S[i,j,k,l] 
SB/: SB[j., k_, L] ** S[i.] := 

S[k, 1] Eta[i, j] + S[l,j] Eta[i, k] + S[j, k] Eta[i, 1] + 
S[i, j, k, 1] 

SB/: SB[i.]**S[j_,k.,L,m_] := Eta[i,j] SB[k,l,m] - 

Eta[i,k] SB[l,m,j] + Eta[i,l] SB[j,k,m] - 

Eta[i,m] SB[j,k,l] + SB[i,j,k,l,m] 

SB/: SB[L]**S[j_,k.,L,m_,n_] := -Eta[i,j] S[k,l,m,n] - 

Eta[i,k] S[l,m,n,j] - Eta[i,l] S[m,n,j,k] - 

Eta[i,m] S[n,j,k,l] - Eta[i,n] S[j,k,l,m] + S[i,j,k,l,ni,n] 

SB/: SB[j.,k_,L,m_,n.]**S[L] := -Eta[i,j] S[k,l,m,n] - 

Eta[i,k] S[l,ni,n,j] - Eta[i,l] S[m,n,j,k] - 

Eta[i,m] S[n,j,k,l] - Eta[i,n] S[j,k,l,m] - S[i,j,k,l,m,n] 

SB/: SB[i.]**S[j_,k.,L,m_,n_,p.] := - SB[i,j,k,l,m,n,p] + 

(Epsiloii[j,k,l,m,n,p,$l[epi], $2[epi], $3[epi], $4[epi]]* 

SB[i,$l[epi], $2[epi], $3[epi], $4[epi++]])/24 

SB/: SB[j.,k_,L,m_,n.,p_,q.]**S[L] := -S[i,j,k,l,ni,n,p,q]- 

(Epsilon[j,k,l,m,n,p,q,$l[epi], $2[epi], $3[epi]]* 

S[i,$l[epi], $2[epi], $3[cpi++]])/6 

SB/: SB[L]**S[i_,k.,L,m_,n_,p.,q_] := -S[i„i,k,l,m,n,p,q]+ 

(Epsilon[j,k,l,m,n,p,q,$l[epi], $2[epi], $3[epi]]* 

S[i,$l[epi], $2[epi], $3[epi++]])/6 

SB/: SB[i.]**S[j_,k.,L,m_,n_,o.,p.,q_] := 

- SB[i,j,k,l,m,n,o,p,q]- 

(Epsiloii[j,k,l,m,n,o,p,q,$l[epi], $2[epi]]* 
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SB[i,$l[epi], $2[epi++]])/2 
SB/: SB[j.,k_,L,m_,n.,o_,p.,q_,r_]**S[L] := 
Epsiloii[i,j,k,l,m,n,o,p,q,r] - 
Epsiloii[j,k,l,m,ii,o,p,q,r,$l[epi]]*S[i,$l[epi++]] 
SB/: SB[L]**S[,i_,k.,L,m_,n_,o.,p.,q_,r_] := 
Epsilon[i,j,k,l,m,n,o,p,q,r] + 
Epsilon[j,k,l,m,n,o,p,q,r,$l[epi]]*S[i,$l[epi++]] 
S/: S[L,j_]**S[k_,L] := -Eta[k,l] S[i,j] - 
Expand[ Par[S[i,j]**S[k]]**SB[l]] 
S/: S[a_,b.,c.,d_]**S[k_,L] := 
-Eta[k,l] S[a,b,c,d] - 
Expand[ Par[S[a,b,c,d]**S[k]]**SB[l]] 
S/: S[a_,b.,c.,d_,e_,L]**S[k_,L] := 
-Eta[k,l] S[a,b,c,d,e,f] - 
Expand[ Par[S[a,b,c,d,e,f]**S[k]]**SB[l]] 
S/: S[a_,b.,c.,d_,e_,L,g.,h_]**S[k_,L] := 
-Eta[k,l] S[a,b,c,d,e,f,g,h] - 
Expand[ Par[S[a,b,c,d,e,f,g,h]**S[k]]**SB[l]] 
S/: S[a_,b.,c.]**S[k.,L] := 
Eta[k,l] S[a,b,c] + 
Expand[ Par[S[a,b,c]**SB[k]]**S[l]] 
S/: S[a_,b.,c.,d_,e_]**S[k_,L] := 
Eta[k,l] S[a,b,c,d,e] + 
Expand[ Par[S[a,b,c,d,e]**SB[k]]**S[l]] 
S/: S[a_,b.,c.,d-,e-,L,g.]**S[k.,L] := 
Eta[k,l] S[a,b,c,d,e,f,g] + 
Expand[ Par[S[a,b,c,d,e,f,g]**SB[k]]**S[l]] 
S/: S[a_,b.,c.,d-,e-,L,g-,h_,m_]**S[k_,L] := 
Eta[k,l] S[a,b,c,d,e,f,g,h,m] + 
Expand[Par[S[a,b,c,d,e,f,g,h,m]**SB[k]]**S[l]] 
S/: S[i.,j_]**S[k_,L.] := -Eta[i,j] S[k,l] - 
Expand[ S[i]**Par[SBp]**S[k,l]]] 
S/: S[i.,j_]**SB[k.,L_] := Eta[i,j] SB[k,l] + 
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Expand[ SB[i]**Par[Sp]**SB[k,l]]] 

SB/: SB[i.,j..]**S[k.,L] := - Eta[k,l] SB[i,j] - 

Expand[Par[SB[i,j]**S[k]]**SB[l]] 

S/: S[i.,j_,k.]**S[l_,ni.,n__] := 

- Expand[S[i]**Par[SB[j]**Par[S[k]**S[l,ni,n]]]] - 
Eta[i,j] S[k]**S[l,m,n] + Eta[i,k] S[j]**S[l,ni,n]- 
Eta[j,k] S[i]**S[l,m,n] 

S/: S[L,j_,k.]**SB[L,m_,n..] := 

- Expand[S[i]**Par[SB[j]** 
Par[S[k]**SB[l,m,n]]]] - 

Eta[i,j] S[k]**SB[l,m,ii] + Eta[i,k] S[j]**SB[l,m,n] - 

Eta[j,k] S[i]**SB[l,ni,n] 

S/: S[i.,j_,k..]**S[L,ni.,n_] := 

- Expand[Par [Par [S[i,j ,k] **S[1]] **SB [m]] ** 
S[n]] - 

Eta[l,m] S[i,j,k]**S[n] + Eta[l,n] S[i,j,k]**S[m] - 

Eta[m,n] S[i,j,k]**S[l] 

S/: S[a_,b.,c..]**SB[L,j_,k.] := 

- Expand[Par[Par[S[a,b,c]**SB[j]]**S[k]]** 
SB[i]] - 

Eta[j,k] S[a,b,c]**SB[i] + Eta[i,j] S[a,b,c]**SB[k] - 

Eta[i,k] S[a,b,c]**SB[j] 

S/: SB[a.,b_,c__]**S[i.,j_,k.] := 

- Expand[Par[Par[SB[a,b,c]**S[i]]**SB[j]]**S[k]] - 
Eta[j,k] SB[a,b,c]**S[i] - Eta[i,j] SB[a,b,c]**S[k] + 
Eta[i,k] SB[a,b,c]**S[j] 

SB/: SB[L,j.,k_]**S[L,m_,n__] := 

- Expand[SB[j]**Par[S[k]** 
Par[SB[i]**S[l,m,n]]]] + 

Eta[i,j] SB[k]**S[l,m,n] - Eta[i,k] SB[j]**S[l,m,n] - 
Eta[j,k] SB[i]**S[l,ni,n] 
S/: S[a_,b.,c.,d_]**S[i.,j_,k.,L.] := 
Expand[S[b,c,d]**Par[SB[a]**S[i,j,k,l]]] + 
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Eta[a,b] Par[S[c,d]**S[i,j,k,l]]- 
Eta[a,c] Par[S[b,d]**S[i,j,k,l]] + 
Eta[a,d] Par[S[b,c]**S[i,j,k,l]] 
S/: S[a_,b.,c.,d-]**SB[i_,j.,k_,l__] := 

- Expand[SB[a]**Par[S[b,c,d]**SB[i,j,k,l]]] + 
Eta[a,b] Par[S[c,d]**SB[i,j,k,l]]- 

Eta[a,c] Par[S[b,d]**SB[i,j,k,l]] + 
Eta[a,d] Par[S[b,c]**SB[i,j,k,l]] 
S/: S[a_,b.,c.,d-]**S[i-,j-,k.,L] := 

- Expand[Par[S[a,b,c,d]**S[i]]**SB[j,k,l]]- 
Eta[i,j] Par[S[a,b,c,d]**S[k,l]]+ 
Eta[i,k] Par[S[a,b,c,d]**Sp,l]]- 

Eta[i,l] Par[S[a,b,c,d]**Sp,k]] 

S/: S[a_,b.,c.,d-,e_,f_]**S[i.,j_,k_,L] := 

- Expand[Par[S[a,b,c,d,e,f]**S[i]]**SB[j,k,l]]- 
Eta[i,j] Par[S[a,b,c,d,e,f]**S[k,l]] + 
Eta[i,k] Par[S[a,b,c,d,e,f]**S[j,l]]- 

Eta[i,l] Par[S[a,b,c,d,e,f]**S[j,k]] 

S/: S[a_,b.,c.,d-,e-,L,g-,h_]**S[L,j_,k.,L] := 

- Expand[Par[S[a,b,c,d,e,f,g,h]**S[i]]**SB[j,k,l]]- 
Eta[i,j] Par[S[a,b,c,d,e,f,g,h]**S[k,l]]+ 
Eta[i,k] Par[S[a,b,c,d,e,f,g,h]**S[j,l]]- 

Eta[i,l] Par[S[a,b,c,d,e,f,g,h]**S[j,k]] 
S/: S[a_,b.,c.,d_,e_]**S[i.,j_,k.,L] := 

- Expand [ 

Par[S[a,b,c,d,e]**SB[i]]**S[j,k,l]] + 
Eta[i,j] Par[S[a,b,c,d,c]**S[k,l]]- 
Eta[i,k] Par[S[a,b,c,d,c]**S[j,l]] + 
Eta[i,l] Par[S[a,b,c,d,e]**S[j,k]] 
S/: S[a_,b.,c.,d_,e_,L,g.]**S[L,j.,k_,L] := 

- Expand [ 

Par[S[a,b,c,d,e,f,g]**SB[i]]**Sp,k,l]]+ 
Eta[i,j] Par[S[a,b,c,d,e,f,g]**S[k,l]]- 
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Eta[i,k] Par[S[a,b,c,d,e,f,g]**S[j,l]]+ 

Eta[i,l] Par[S[a,b,c,d,e,f,g]**Sp,k]] 

S/: S[a_,b.,c.,d_,e_,f_,g.,h_,m_]**S[i.,j_,k.,L] := 

- Expand [ 

Par[S[a,b,c,d,e,f,g,h,m]**SB[i]]**S[j,k,l]] + 
Eta[i,j] Par[S[a,b,c,(Lc,f,g,h,ni]**S[k,l]]- 
Eta[i,k] Par[S[a,b,c,d,c,f,g,h,ni]**S[j,l]] + 
Eta[i,l] Par[S[a,b,c,d,e,f,g,h,m]**S[j,k]] 
S/: SB[a.,b_,c_,d..]**S[i_,j.,k_,l_] := 

- Expand[Par[SB[a,b,c,d]**S[i]]**SB[j,k,l]]- 
Eta[i,j] Par[SB[a,b,c,d]**S[k,l]] + 
Eta[i,k] Par[SB[a,b,c,d]**Sp,l]]- 

Eta[i,l] Par[SB [a,b,c,d]**Sp,k]] 
S/: S[a_,b.,c.,d_,e_]**SB[i_,j.,k_,l_,m_] := 
Expand[S[a]**Par[S[b,c,d,e]**SB[i,j,k,l,m]]]- 
Eta[a,b] Par[S[c,d,e]**SB[i,j,k,l,m]] + 
Eta[a,c] Par[S[b,d,c]**SB[i,j,k,l,m]] - 
Eta[a,d] Par[S[b,c,c]**SB[i,j,k,l,m]] + 
Eta[a,c] Par[S[b,c,d]**SB[i,j,k,l,in]] 
SB/: SB[a_,b.,c.,d.,e_]**S[L,j_,k.,L,ni.] := 
Expand[S[b,c,d,e]**Par[SB[a]**S[i,j,k,l,m]]] + 
Eta[a,b] Par[SB[c,d,e]**S[i,j,k,l,ni]] - 
Eta[a,c] Par[SB[b,d,e]**S[i,j,k,l,ni]] + 
Eta[a,d] Par[SB[b,c,e]**S[i,j,k,l,ni]] - 
Eta[a,e] Par[SB[b,c,d]**S[i,j,k,l,ni]] 
Unprotect [Plus] 

x_?NumbcrQ Epsilon[a__,b__] S_[b__] + 
y_?NumbcrQ Epsilon[a__,c__] S_[c__] := 
(x+y) Epsilon[a,b]**S[b] 

x_?NumberQ Epsilon[l ,m_,n__,o__] S_[m_,o__] 

y_?NuniberQ Epsiloii[l ,m_,n__,p__] S_[ni_,p__] 

(x+y) Epsilon[l,ni,n,o]**S[m,o] 
x_?NuniberQ Epsilon[a_.,b__] S_[i_,b._] + 



y_?NuniberQ Epsilon[a__,c__] S_[i_,c__] := 
(x+y) Epsilon[a,b] S[i,b] 

x_?NuniberQ Epsilon[l ,m_,ii__,o__] S_p_,ni_,o__] + 

y_?NuniberQ Epsilon[l ,m_,n__,p__] S_p_,ni_,p__] := 

(x+y) Epsilon[l,m,n,o] S[j,m,o] 

x_?NumbcrQ Epsilon[l ,m_,n__,o__] S_[m_,j_,o__] + 

y_?NumbcrQ Epsilon[l ,m_,n__,p__] S_[m_,j_,p__] := 

(x+y) Epsilon[l,m,n,o] S[m,j,o] 

Protect [Plus] 

(* EndPackageO *) 
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